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Porous  carbon  materials  are  typically  used  in  both  the  substrate  (typically  carbon  paper)  and  the  elec¬ 
trocatalyst  supports  (often  platinised  carbon)  within  proton  exchange  membrane  fuel  cells.  Gravimetric 
nitrogen  adsorption  has  been  studied  at  a  carbon  paper  substrate,  two  different  Pt-loaded  carbon  paper 
electrodes  and  three  particulate  carbon  blacks.  N2  BET  surface  areas  and  surface  fractal  dimensions  were 
determined  using  the  fractal  BET  and  Frenkel-Halsey-Hill  models  for  all  but  one  of  the  materials  studied. 
The  fractal  dimensions  of  the  carbon  blacks  obtained  from  gas  adsorption  were  compared  with  those 
obtained  independently  by  small  angle  X-ray  scattering  and  showed  good  agreement.  Density  functional 
theory  was  used  to  characterise  one  of  the  carbon  blacks,  as  the  standard  BET  model  was  not  applicable. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  varieties  of  carbon  [1]  include  activated  carbon,  carbon 
black,  carbon  fibres  and  carbon  cloth,  with  clear  definitions  of  these 
and  all  the  forms  of  carbon  as  a  solid  described  by  Fitzer  et  al.  [2], 
In  most  catalytic  applications,  including  proton  exchange  mem¬ 
brane  (PEM)  fuel  cells  [3],  carbon-containing  composites  are  the 
substrates  under  consideration.  Carbon  paper  is  often  chosen  as 
the  substrate  in  PEM  fuel  cells  and  is  usually  prepared  from  the 
random  weaving  of  carbon  fibres,  with  the  extent  and  nature  of 
the  surface  area  and  porosity  being  dependent  upon  the  precursor 
material  and  its  manufacturing  history. 

Carbon  blacks  are  often  used  as  the  support  material  for  pre¬ 
cious  metal  electrocatalyst  coatings.  Carbon  blacks  play  a  major 
role  in  structural  components,  conductive  supports,  and  electrocat¬ 
alysts  in  batteries  and  fuel  cells.  Such  materials  possess  a  relatively 
high  surface  area  which  is  derived  from  the  fine  particle  size  of  the 
material,  rather  than  from  porosity.  Optimising  the  structural  prop¬ 
erties  of  carbon  blacks  is  an  important  area  of  research  due  to  their 
importance  in  fuel  cell  applications. 
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The  characterisation  of  carbon  blacks  has  previously  been 
determined  using  a  comparison  of  nitrogen  and  carbon  diox¬ 
ide  adsorption  at  77  and  298  K,  respectively  [4],  atomic  force 
microscopy  [5],  Raman  and  neutron  scattering  [6],  X-ray  diffraction 
[7],  small  angle  X-ray  scattering  (SAXS)  [8],  transmission  electron 
microscopy  and  small  angle  neutron  scattering  [9],  Xu  et  al.  [10] 
used  various  gas  adsorbates  to  investigate  the  structural  proper¬ 
ties  of  different  grades  of  carbon  blacks.  Analysis  of  the  data  using 
the  work  of  Pfeifer  and  Avnir  [11  ]  and  Halsey  [12]  determined  the 
fractal  dimension  of  the  materials. 

Carbon  fibre  materials  tend  to  have  carbon  content  greater  than 
90%.  Since  the  final  carbon  fibres  contain  almost  100%  carbon,  any 
fibrous  material  with  a  carbon  backbone  can  potentially  be  used  as 
a  precursor  which  would  yield  a  carbonaceous  residue.  Rayon  and 
polyacrylonitrile  are  precursors  for  many  of  the  commercial  fibres 
that  are  produced,  although  other  precursors  such  as  pitch  (a  by¬ 
product  of  the  petroleum  or  coal-coking  industry),  phenolic  resins 
and  polyacetylenes  have  also  been  used  [13].  Gas  adsorption  is  a 
common  technique  in  the  study  of  structural  properties  of  carbon 
fibres;  a  review  by  Do  and  Do  [14]  has  considered  the  adsorp¬ 
tion  of  supercritical  fluids  (Kr,  Ar,  N  and  CH4)  on  different  porous 
and  non-porous  carbons.  Studies  using  Kr  [15-17]  and  N2  [18,19] 
have  also  been  performed.  Techniques  such  as  scanning  electron 
microscopy  and  atomic  force  microscopy  [20],  X-ray  diffraction  and 
Raman  spectroscopy  [21],  scanning  tunnelling  microscopy  [22,23] 
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Nomenclature 

C  BET  constant 

d  fractal  dimension 

Ej  energy  of  adsorption 

El  liqueflcation  energy 

n  amount  of  gas  adsorbed  (mol) 

nm  molecular  capacity 

Na  Avagadro  constant  (6.023  x  1023  mol-1 ) 

P  applied  pressure  (bar) 

Po  saturation  pressure  (bar) 

R  molar  gas  constant  (8.314J  mol-1  K-1 ) 

T  temperature  (K) 

V  volume  of  gas  adsorbed  at  equilibrium  pressure 

(m3) 

Vm  volume  of  gas  in  monolayer  (m3 ) 

Greek  letter 

a  molecular  cross-sectional  area  (nm2 ) 


and  SAXS  [24]  have  been  used  to  investigate  the  structure  and  phys¬ 
ical  properties  of  carbon  fibres. 

The  work  presented  here  aims  to  interpret  isothermal  data, 
using  a  selection  of  established  models,  for  a  variety  of  impor¬ 
tant  carbon  materials  that  are  typically  and  routinely  used  in  PEM 
fuel  cells.  The  materials  include  a  carbon  paper  substrate  for  elec¬ 
trocatalysts  and  also  two  Pt-loaded  electrodes.  In  addition,  three 
carbon  blacks,  which  are  used  as  the  electrocatalyst  support,  are 
examined.  The  results  aim  to  give  a  quantitative  and  accurate 
description  of  the  structural  properties  of  these  materials.  Both 
classical  and  recently  developed  adsorption  models  are  considered 
and  the  paper  represents  an  ambitious  attempt  to  fit  these  models 
to  standard  isotherms  of  practical  PEM  fuel  cell  carbons.  Features 
of  this  paper  include:  (a)  a  quantitative  consideration  of  several 
adsorption  models,  (b)  classical  adsorption  models  are  comple¬ 
mented  by  modern,  (e.g.,  modified  fractal)  ones,  (c)  adsorption  data 
on  carbon  materials  in  the  paper  have  not  previously  been  reported, 
(d)  not  all  of  the  materials  can  be  fitted  to  simple  adsorption  models 
and  (e)  problems  in  applying  adsorption  models  to  certain  carbons 
are  highlighted. 

2.  Theory 

2.1.  Determination  of  surface  area— the  BET  method 


made  that  E,-  is  not  dependent  on  surface  coverage,  0,-.  This  implies 
both  a  uniform  array  of  surface  sites  and  the  absence  of  lateral  inter¬ 
actions  between  the  adsorbed  molecules.  A  further  assumption  is 
made  to  the  model,  which  is  that  the  multilayer  has  an  infinite 
thickness  at  P/P0  =  1.  This  assumption  allows  the  simplification  of 
the  summation  of  the  amounts  adsorbed  in  all  layers  and  allowed 
the  production  of  the  BET  equation: 

_ P/Pa  _  JL  +  —  (2) 

n(l  -  P/P0)  nmC  +  nmc  P0 

where  nm  is  the  monolayer  capacity  and  C  is  an  empirical  constant. 

According  to  the  BET  theory,  the  constant  C  is  related  exponen¬ 
tially  to  Ej  by  the  simplified  equation: 


C  exp 


(3) 


where  the  quantity  E\  —  El  is  termed  the  net  molar  energy 
of  adsorption.  The  value  of  C  is  an  indication  of  the 
adsorbent-adsorbate  interactions,  with  high  values  (0100) 
indicative  of  high  interaction  strengths,  whilst  medium  (<80)  and 
low  values  (<20)  represent  medium  and  low  interaction  strengths. 

The  validity  of  the  original  BET  equation  is  always  confined  to  a 
limited  part  of  the  isotherm,  which  is  seldom  above  P/P0  =  0.3.  The 
BET  model  does  not  take  into  account  lateral  interactions  between 
the  adsorbed  molecules  in  the  first  layer,  and  assumes  that  all 
higher  layers  have  liquid-like  properties.  Another  prerequisite  of 
the  model  is  that  the  surface  sites  are  all  equivalent  (i.e.  a  uni¬ 
form  surface)  and  that  the  surface  is  flat.  Although  the  BET  method 
includes  assumptions  and  limitations,  it  is  widely  used  for  the  eval¬ 
uation  of  the  surface  area  from  physisorption  isotherm  data. 


2.1.2.  Determination  of  surface  area 

Two  stages  are  involved  in  the  evaluation  of  the  surface  area. 
Firstly,  it  is  necessary  to  construct  the  BET  plot  and  to  derive  a 
value  of  nm.  The  next  stage  is  the  calculation  of  the  specific  surface 
area,  a(BET),  which  requires  knowledge  of  the  average  area,  a,  occu¬ 
pied  by  each  molecule  in  the  completed  monolayer.  The  value  of  nm 
and  C  can  be  solved  by  plotting  (P/P0)/[n(l  -  P/Pq  )]  against  P/Pq.  The 
usual  range  of  linearity  of  the  BET  plot  is  0.1  <P/Pq  <  0.3.  This  range 
varies,  depending  on  the  material  investigated.  Sufficient  experi¬ 
mental  points  on  the  adsorption  isotherm  must  be  determined  to 
select  the  region  providing  the  best  linear  fit.  The  second  stage  in 
the  application  of  the  BET  method  is  the  calculation  of  the  specific 
surface  area,  a(BET),  from  the  monolayer  capacity  nm: 


2.1.1.  The  BET  model 

The  most  popular  technique  for  the  determination  of  the  sur¬ 
face  area  over  a  wide  range  of  porous  materials  (the  BET  model)  is 
the  method  proposed  by  Brunauer  et  al.  [25].  According  to  the  BET 
model,  the  molecules  in  one  layer  can  act  as  possible  sites  for  the 
adsorption  of  molecules  in  the  next  layer.  When  the  rate  of  conden¬ 
sation  (adsorption)  is  equal  to  the  rate  of  evaporation  (desorption), 
an  equilibrium  pressure,  P,  is  achieved  at  a  given  temperature,  T. 
For  adsorption  in  each  layer,  the  following  expression  applies: 

aiPei.^bA  ex p(-J|)  (1) 

where  0;_ i  and  0;  are,  respectively,  the  fractions  of  surface  covered 
by  the  i  — 1  and  i  layers,  a*  and  b,  are  adsorption  and  desorption 
constants  and  E,  is  the  energy  of  adsorption  for  the  ith  layer.  Each 
layer  has  a  different  set  of  values  of  a„  b,  and  E,-  but  it  is  assumed 
that  for  all  layers  after  the  first,  these  parameters  remain  constant 
and  also  that  E,  =  EL,  the  liquefaction  energy.  Also,  the  assumption  is 


a(BET)  =  nmNA  a  (4) 

where  NA  is  the  Avagadro  constant  (6.023  x  1023  mol-1 )  and  a  is 
the  average  area  occupied  by  each  adsorbed  molecule  in  the  com¬ 
pleted  monolayer.  Varied  values  of  a  are  found  in  the  literature; 
the  value  of  a  is,  to  some  extent,  dependent  on  the  adsorbent- 
adsorptive  system  and  temperature  and  not  just  on  the  adsorbate. 
If  the  adsorbent-adsorbate  interactions  are  strong  and  the  tem¬ 
perature  is  low,  the  monolayer  is  likely  to  be  localised  and,  under 
these  conditions,  the  monolayer  structure  tends  to  be  controlled  by 
surface  chemistry. 

The  most  widely  used  adsorptive  for  surface  area  determina¬ 
tion  is  nitrogen  (at  77  K)  as  it  is  considered  to  be  the  most  suitable 
adsorptive  for  the  determination  of  the  surface  area  of  nonporous, 
macroporous,  or  mesoporous  solids.  It  is  assumed  that  the  BET 
nitrogen  monolayer  is  close-packed,  giving  a  <t(N2)  =  0.162  nm2  at 
77  K  [26]. 
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2.2.  Fractal  analysis 

The  concept  of  fractality  and  the  definition  of  the  fractal  dimen¬ 
sion  are  related  to  the  nature  of  the  measurements  of  the  surface 
area  and  other  geometrical  parameters  of  porous  solids.  The  fractal 
approach  provides  a  method  of  evaluating  the  absolute  area  of  a 
finely  divided  or  porous  solid. 

The  fractal  dimension  is  a  structural  parameter  that  quantifies 
the  scaling  invariance  of  self-similar  systems  [27],  If  an  object  can 
be  decomposed  into  X  similar  parts  of  size  n  times  smaller  than 
the  whole,  the  fractal  dimension  is  defined  as  a  scaling  exponent 
in  the  relation  between  the  fragmentation  parameter  X  and  the 
contraction  ratio,  n: 

A  =  nd  (5) 

The  fractal  dimension  of  self-similar  objects  can  be  determined 
by  different  methods.  In  the  method  of  “box  counting”,  the  fractal 
is  covered  by  equal  boxes  of  size  a.  The  number  of  boxes,  N  depends 
on  the  box  size  a,  and  the  volume  V occupied  by  the  boxes  and  scales 
as: 

V(a)  a  a3~d  (6) 

The  method  of  molecular  tiling  was  introduced  by  Pfeifer  and 
Avnir  [11]  and  involves  a  comparison  of  the  monolayer  capacities 
of  different  adsorbates.  The  molecules  of  different  sizes  are  used  as 
‘gauges’  to  determine  the  surface  area  of  an  adsorbent.  In  a  stan¬ 
dard  adsorption  experiment,  the  surface  area,  S,  is  estimated  as  the 
product  of  the  number  of  molecules  in  the  monolayer  capacity  nm 
and  the  effective  area  occupied  by  one  molecule,  a: 

S  =  nmcr  (7) 

According  to  the  method  of  tiling,  the  number  of ‘gauges’  needed 
to  cover  a  fractal  surface  scale  with  the  linear  size  of  the  gauges, 
a  =  cr1/2  and: 


where  the  subscript  ‘0’  denotes  the  adsorbate  taken  as  a  refer¬ 
ence.  The  fractal  dimension,  d,  can  be  determined  from  a  set  of 
adsorption  measurements  using  different  adsorbates.  The  mono- 
layer  capacity  is  calculated  from  the  gas  adsorption  isotherm  using 
the  BET  method,  as  mentioned  earlier.  A  problem  of  using  the 
method  of  molecular  tiling  is  the  limited  range  of  length-scales 
available  for  fractal  analysis.  The  method  can  reveal  information 
about  the  molecular  scale  roughness  of  the  adsorbent  surface  but 
it  is  difficult  to  make  conclusions  about  the  fractality,  which  implies 
a  hierarchical  structure  over  a  wide  range  of  scales. 

Pfeifer  et  al.  [28,29]  adopted  the  Frenkel-Halsey-Hill  (FHH)  [12] 
equation  for  multilayer  adsorption  of  fractal  surfaces.  The  method 
is  based  on  an  expression  for  the  surface  fractal  dimension  from  an 
analysis  of  multilayer  adsorption  on  a  fractal  surface,  such  that: 

ln(iD-c+slnK£)]  (9) 

where  V  is  the  volume  of  gas  adsorbed  at  equilibrium  pressure  P, 
Vm  is  the  volume  of  gas  in  a  monolayer,  and  P0  is  the  saturation 
pressure.  The  constant  C  is  a  pre-exponential  factor,  and  S  is  a  power 
law  exponent  dependent  on  d,  the  surface  fractal  dimension  and 
the  mechanism  of  adsorption.  There  are  two  limiting  cases  for  the 
model.  At  the  lower  end  of  the  isotherm,  which  represents  the  early 
stages  of  multilayer  build-up,  the  film/gas  interface  is  controlled 
by  attractive  van  der  Waal’s  forces  between  the  gas  and  the  solid. 
This  tends  to  allow  the  film/gas  interface  to  replicate  the  surface 


roughness.  The  value  of  the  constant,  S  is  given  by: 


S  = 


d-3 

3 


(10) 


When  there  is  a  high  degree  of  coverage,  however,  the  interface 
is  determined  by  the  liquid/gas  surface  tension  which  causes  the 
interface  to  move  further  away  from  the  surface  and  reduces  the 
surface  area.  In  this  case,  S  is  given  by: 


(11) 


The  surface  fractal  dimension,  d  may  physically  take  values  in 
the  range  2  <  d  <  3,  with  a  value  of  2  indicating  a  perfectly  smooth 
surface,  whilst  a  value  of  3  is  indicative  of  a  highly  disordered 
surface.  An  alternative  analysis  method  to  derive  surface  fractal 
dimensions  from  gas  adsorption  data  is  to  use  the  fractal  version 
of  the  BET  theory.  The  concept  is  based  upon  the  fact  that,  for  a 
fractal  surface,  the  area  available  for  adsorption  in  the  ith  layer  of 
adsorptive  decreases  by  the  factor/]  and  is  given  by: 

Ai 

=  At  = 

where  a  =  3-d.  Mahnke  and  Mogel  [30]  suggested  an  alternative 
expression  to  the  fractal  BET  equation  initially  proposed  by  Fripiat 
et  al.  [31]: 

LSI  w  =  [i  -(P^ctP/Po)] 


(12) 


-a  log[l-(P/P0)]  (13) 

where  C  is  the  BET  constant.  Mahnke  and  Mogel  [30]  suggested 
that  the  following  procedure  should  be  performed  to  obtain  a  value 
for  a:  (a)  estimate  C  and  l/m  from  the  BET  [25]  equation  then  (b) 
evaluate  the  plot  of: 

«•  -«-('« 
(M) 


The  plot  should  be  linear  and  the  slope  provides  an  estimate  for 
a  =  3-d. 

It  has  previously  been  shown  [32]  that  adsorption  of  different 
adsorbates  on  heterogeneous  surfaces  may  vary  due  to  interac¬ 
tion  strength  differences  between  the  adsorbate  and  adsorbent.  In 
a  similar  way,  if  a  surface  is  composed  of  multiple  components, 
these  interaction  strength  differences  may  also  occur.  This  can  be 
investigated  by  using  a  two-component  BET  model  [32]: 

V  P/Po  \  flCi  (1-0)C2  1 

i  -  m)3_d  L1 + (ci  -  ’W)  i + (Q  -  i)(p/p0)J 

(15) 


The  overall  isotherm  is  thus  a  composite  of  two  BET-type  equa¬ 
tions,  each  one  corresponding  to  the  two  types  of  surface  patches. 
Cj  and  C2  are  the  BET  constants  for  surface  patches  of  types  1  and 
2,  respectively,  d  is  the  surface  fractal  dimension  and  0  is  the  frac¬ 
tion  of  the  surface  occupied  by  patches  of  type  1.  This  approach 
may  highlight  the  differences  in  strength  of  interaction  between 
the  adsorbed  nitrogen  and  the  carbon  blacks  and  Pt  within  the  fuel 
cell  electrodes. 

2.3.  Small  angle  X-ray  scattering 

Small-angle  X-ray  scattering  (SAXS)  is  a  technique  that  allows 
the  determination  of  the  particle  size  distribution  and  the  mea¬ 
surement  of  the  specific  surface  area.  A  beam  of  X-rays  is  scattered 
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by  the  electrons  in  an  irradiated  material.  As  the  distribution  of 
electrons  throughout  a  material  is  heterogeneous,  fluctuations  in 
electron  density  (number  of  electrons  per  unit  volume)  exist  at 
different  zones  within  the  sample.  When  the  dimensions  of  these 
spatial  inhomogeneities  are  similar  to  the  wavelength  of  the  inci¬ 
dent  X-ray  beam,  X  then  most  of  the  scattering  is  observed  at  angles 
greater  than  10°.  However,  if  the  inhomogeneities  in  the  electron 
density  extend  over  the  distances  in  the  range  between  ca.  0.5 
and  400  nm,  then  the  intensity  of  scattered  X-rays  is  appreciable 
at  small  angles.  The  SAXS  technique  can  be  used  to  provide  infor¬ 
mation  on  structures  much  larger  than  the  normal  inter-atomic 
distances  encountered  in  dense  materials. 

Porous  materials  greatly  scatter  X-ray  radiation  near  the  X-ray 
beam,  due  to  the  system  of  volume-distributed  local  inhomo¬ 
geneities  of  electron  density  (submicropores  and  micropores).  Pore 
sizes  and  shapes,  exceeding  the  wavelength  of  X-ray  radiation, 
condition  the  angular  range  of  such  scattering  and  the  scatter¬ 
ing  intensity,  I(q)  is  mainly  determined  by  pore  concentration  and 
electron  density  gradient  at  the  pore-matrix  border. 

The  surface  roughness  on  a  scale  in  the  range  of  the  inverse 
of  q  covered  in  the  SAXS  measurements  can  be  determined  [33], 
This  is  achieved  by  describing  the  material  as  having  surface  fractal 
properties,  and  a  decay  in  the  scattered  intensity  I(q)  is  given  by: 

J(<j)~<rC6-d)  (i6) 


where  d  is  the  surface  fractal  dimension.  Typically,  values  of  2  <  d  <  3 
occur,  with  d  =  2  representing  a  smooth  surface,  whilst  d  =  3  is 
indicative  of  a  highly  disordered  surface.  The  surface  fractal  dimen¬ 
sion  obtained  from  SAXS  can  also  be  used  as  the  basis  to  obtain, 
independently  of  gas  adsorption,  an  estimate  of  the  specific  sur¬ 
face  area,  as  perceived  by  a  given  adsorbate,  of  a  porous  solid  with 
a  fractally  rough  internal  surface.  If  the  porous  solid  consists  of  a 
close-packing  of  spherical  particles  of  radius  R,  which  each  possess 
a  fractally  rough  surface  of  dimension  d,  then  the  specific  surface 
area  of  the  solid  (As)  as  perceived  by  a  particular  adsorbate  species 
is  given  by: 


„  3  r2(R/r)d 

A^=-^r 


(17) 


where  r2  is  the  cross-sectional  area  of  a  close-packed  adsorbate 
molecule  and  ps  is  the  skeletal  mass  density  of  the  solid.  For  porous 
solids  consisting  of  packed  particles,  the  characteristic  particle 
size  (R)  can  be  obtained  using  electron  microscopy  [33],  mercury 
porosimetry  [34]  or  SAXS  [35], 


3.  Experimental  details 

The  carbon  materials  examined  in  this  work  include  those  con¬ 
ventionally  used  within  a  membrane  electrode  assembly  (MEA)  in  a 
PEM  fuel  cell.  Three  types  of  carbon  paper  were  examined,  namely 
(i)  an  unloaded  (denoted  CU1 ),  (ii)  a  low  Pt-loading  (0.39  mg  cm-2, 
denoted  CL1)  and  (iii)  a  high  Pt-loading  (4.7  mg  cm-2,  denoted 
CHI ).  The  PEM  fuel  cell  electrodes  consisted  of  a  carbon  paper  sub¬ 
strate  and  a  mixture  of  activated  carbon  and  Pt  as  the  support  and 
catalyst,  respectively.  Electrode  materials  were  supplied  by  John¬ 
son  Matthey  Technology  Centre.  Both  CHI  and  CL1  were  supported 
by  TGP-H-060  Toray  paper.  The  catalyst  layer  of  CL1  was  made 
up  of  40  wt%  Pt  supported  on  Vulcan  XC72R  containing  75wt% 
Nafion®  1100  EW  ionomer  (based  on  carbon  weight)  at  a  loading  of 
0.39  mg  Pt  cm-2.  The  catalyst  layer  of  CHI  was  made  of  Pt  black  con¬ 
taining  6  wt%  Teflon  at  4.4  mg  Pt  cm-2 .  Both  samples  also  had  a  top 
coat  of  Nafion®  ion  exchange  resin  at  a  loading  of  0.5  mg  Pt  cm-2. 

Three  types  of  particulate  carbon  used  as  a  Pt  support  material, 
were  also  examined.  These  included  (1)  a  conductive  carbon,  (2) 


Vulcan  XC72R  and  (3)  an  acetylene  black  (denoted  CPI,  CP2  and 
CP3,  respectively).  The  structural  and  physical  characteristics  of  the 
materials  were  not  known  prior  to  analysis. 

3.1.  Experimental  methods 

3.1.1.  Scanning  electron  microscopy 

Scanning  electron  micrographs  were  obtained  for  the  carbon 
paper  samples  (CU1,  CL1  and  CHI)  using  a  JEOL JSM-6310  micro¬ 
scope.  The  samples  were  cut,  using  a  clean  steel  scalpel  blade,  into 
0.5  cm2  squares  and  fixed  to  an  aluminium  sample  dish  using  a 
conductive  carbon  tape.  The  imaging  was  conducted  at  an  operat¬ 
ing  voltage  of  15  kV.  Images  were  acquired  and  stored  directly  to  a 
PC  using  the  image  capture  software  Semafore  version  4.  Multiple 
samples  of  the  same  batch  were  examined  to  investigate  structural 
heterogeneity. 

3.1.2.  Nitrogen  adsorption 

Nitrogen  adsorption  experiments  were  performed  at  77  K  using 
a  Micromeritics  accelerated  surface  area  and  porosimetry  (ASAP) 
2010  apparatus.  Carbon  paper  samples  were  cut  with  a  clean  scalpel 
blade  into  50  strips,  each  having  dimensions  of  0.7  cm  x  2  cm.  The 
number  of  strips  was  kept  constant  in  all  experiments.  The  sample 
(collection  of  50  strips)  was  then  placed  into  a  pre-weighed  round- 
bottomed  glass  tube  and  a  glass  spacer  was  placed  gently  into  the 
sample  tube  to  reduce  the  free-space  volume.  A  sealing  frit  was 
then  fitted  to  the  opening  of  the  tube  so  as  to  seal  the  tube  after 
degassing  prior  to  analysis.  The  sample  tube  and  its  contents  were 
then  loaded  into  the  degassing  port  of  the  apparatus.  A  heating 
jacket  was  applied  to  the  tube  and  the  sample  was  heated,  under 
vacuum,  to  100  °C  for  24  h. 

The  purpose  of  the  thermal  pre-treatment  for  each  sample  was 
to  drive  off  any  physisorbed  water  on  the  sample  whilst  leaving 
the  morphology  of  the  sample  unchanged.  Once  the  preparation 
was  completed,  the  heating  jacket  was  removed  and  the  sample 
was  allowed  to  cool  down  to  room  temperature  (295  K).  The  sam¬ 
ple  tube  and  its  contents  were  then  re-weighed  to  obtain  the  dry 
weight  of  the  sample,  which  was  entered  into  the  necessary  section 
of  the  software  program,  prior  to  transfer  of  the  sample  tube  to  the 
analysis  port  for  the  automated  analysis  procedure. 

The  sample  tube  was  immersed  in  liquid  nitrogen  at  77  K 
before  the  adsorption  measurements  were  taken  in  the  rela¬ 
tive  pressure  region  of  0.05-1.00  P/P0.  The  software  controlling 
the  automated  apparatus  performed  a  leak-checking  procedure 
and  an  equilibration  time  of  45  s  was  used  for  each  adsorption 
point.  Once  the  experiment  was  completed,  the  data  values  of  the 
adsorption/desorption  isotherms  were  available  for  analysis.  The 
cross-sectional  area  of  nitrogen,  a,  was  taken  to  be  0.162  nm2.  The 
same  procedure  was  used  for  the  carbon  black  materials,  the  nom¬ 
inal  amount  of  each  sample  analysed  being  250  mg. 

3.1.3.  Small  angle  X-ray  scattering 

SAXS  patterns  were  run  on  the  carbon  black  samples  (CP1-CP3) 
using  the  DUBBLE  SAXS  beamline  on  BM26B  at  the  European  Syn¬ 
chrotron  Radiation  Facility  (ESRF).  The  two-dimensional  patterns 
were  corrected  for  detector  efficiency  and  calibrated  against  silver 
behenate  and  pixels  at  the  same  radius  were  averaged  to  produce 
one-dimensional  patterns. 

3.1.4.  Transmission  electron  microscopy 

Transmission  electron  micrographs  of  the  carbon  paper  samples 
(CU1,  CL1  and  CHI )  were  obtained  using  ajeol  1200  microscope. The 
samples  were  embedded  in  a  spurr’s  resin  and  vacuum-treated  to 
enhance  resin  penetration  and  to  ensure  that  the  samples  were 
degassed.  The  samples  were  sectioned  using  a  microtome  and 
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freshly  prepared  glass  knives.  Images  were  acquired  using  a  35  mm 
SLR  camera,  the  negatives  being  digitized  for  analysis. 

4.  Results  and  discussion 

4.1.  Carbon  paper  samples  CU1-CH1 

Fig.  1  shows  a  series  of  scanning  electron  micrograph  (SEM) 
images  of  the  carbon  paper  electrodes  investigated  in  this  work. 
Fig.  1(a)  is  a  micrograph  of  CHI  and  is  comprised  of  a  series  of 
interwoven  and  overlapping  carbon  fibres,  with  typical  diameters 
over  the  range  of  5-15  p,m.  A  micrograph  of  an  individual  fibre  is 
shown  in  Fig.  1(b)  and  illustrates  the  ‘plate-like’  structure  of  the 
fibre  together  with  the  presence  of  grooves  along  the  axial  dimen¬ 
sion.  Fig.  1(a)  shows  that  the  carbon  substrate  does  not  possess 
pores,  as  defined  in  the  conventional  sense  but  open  areas;  ‘holes’ 
exist  between  the  overlapping  fibres.  These  holes  vary  in  size  from 
ca.  10  to  100  pun,  whilst  Fig.  1(b)  implies  that  the  individual  fibres 
have  limited  or  zero  porosity  themselves  (based  on  the  scales  acces¬ 
sible  by  electron  microscopy). 

Fig.  1(c)  and  (d)  are  SEM  images  of  the  carbon  substrate  loaded 
with  0.39  mg  Pt  cm-2  (denoted  CL1).  Fig.  1(c)  shows  how  the  sur¬ 
face  of  CL1  consists  of  irregular-shaped  ‘packets’  of  deposited  Pt, 
with  large  cracks  seen  between  the  individual  packets.  It  can  be  fur¬ 
ther  seen  in  Fig.  1(c)  that  small  imperfections  in  the  Pt  layer  have 
caused  small  holes  and  larger  agglomerates  randomly  positioned 
within  the  surface.  Fig.  1(c)  and  (d)  illustrates  a  certain  degree  of 
heterogeneity  in  terms  of  surface  deposition  and  coverage.  This 
heterogeneity  is  obtained  due  to  the  larger  overall  volume  of  the 
catalyst  used  in  the  preparation  of  the  electrode.  Fig.  1(d)  shows  an 
individual  fibre  that  is  coated  in  Pt  and  is  at  a  considerable  depth 
within  the  paper  (relative  to  the  overall  thickness  of  the  paper).  This 
sample  is  typically  used  as  an  electrode  in  a  MEA  and  the  purpose  of 


the  Pt  layer  is  to  catalyse  the  oxidation  of  hydrogen  or  the  reduction 
of  oxygen.  If  the  Pt  layer  does  not  cover  the  substrate  effectively,  the 
catalytic  activity  of  the  electrode  decreases;  the  performance  of  the 
electrode  assembly  is  then  lower. 

Fig.  1(e)  shows  a  SEM  micrograph  of  sample  CHI.  The  Pt  loading 
of  CHI  is  over  10  times  greater  than  CL1,  as  its  typical  application  is 
in  membrane  electrode  assemblies  that  require  a  high  current  den¬ 
sity  together  with  a  long  service  lifetime.  Fig.  1(e)  clearly  shows 
that  the  deposition  of  Pt  is  not  uniform,  and  that  areas  with  none 
or  little  Pt  are  present.  This  electrode  has  a  thinner  layer  of  cata¬ 
lyst  (platinum  black)  due  to  the  absence  of  a  carbon  support.  This 
fact  results  in  a  lower  volume  of  material  available  to  coat  the  car¬ 
bon  fibre  substrate.  Fig.  1(f)  shows  the  nodular  arrangement  of  the 
Pt-activated  C  agglomerates,  which  range  from  ca.  10  to  100  nm 
in  size  and  qualitatively  imply  a  higher  surface  area  than  sample 
CL1. 

Fig.  2  shows  a  series  of  transmission  electron  micrographs 
(TEMs)  of  the  carbon  paper  samples  investigated.  Fig.  2(a)  is  a 
TEM  micrograph  of  an  individual  fibre  from  sample  CU1,  and  very 
clearly  illustrates  the  plate-like  formation  of  the  carbon  structure. 
The  TEMs  suggest  a  limited  degree  of  porosity  is  present  with  any 
available  porosity  only  accessible  between  the  plates  themselves. 

TEMs  for  sample  CL1  are  shown  in  Fig.  2(b)  and  (c).  Fig.  2(b) 
shows  two  open,  elliptical  regions  in  the  top-right  and  bottom-left 
of  the  image  and  can  be  attributed  to  two  individual  fibres.  The 
majority  of  the  features  seen  in  Fig.  2(b)  are  the  activated  carbon  and 
Pt  mixture  that  serves  as  the  porous  support  and  electrocatalyst, 
respectively.  Fig.  2(c)  is  a  high-resolution  image  focussing  on  the 
activated  carbon  and  Pt  particles.  The  idealised  structure  for  the 
support  [36]  is  envisaged  as  large  carbon  agglomerates  supporting 
a  highly  dispersed  level  of  a  low  concentration  of  Pt.  The  darkened 
areas  of  Fig.  2(c)  are  the  activated  carbon  support,  whilst  the  lighter 
more  dispersed  and  smaller  parts,  being  the  Pt  catalyst. 
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Fig.  2.  Transmission  electron  micrographs  of  particulate  carbon  samples:  (a)  an  individual  fibre  from  unloaded  sample  CU1;  (b)  and  (c)  low  loaded  (0.39mgcm-2  Pt)  sample 
CL1,  (d)  and  (e)  high  loaded  (4.7  mg  cm-2  Pt)  sample  CHI. 


TEMs  of  sample  CHI  are  shown  in  Fig.  2(d)  and  (e).  Individual 
fibres,  shown  by  the  round  areas  within  the  middle  of  Fig.  2(d)  are 
similar  to  that  seen  in  Fig.  2(b).  The  light  grey  area  is  the  resin  used 
to  set  the  sample,  whilst  the  dark  black  areas  are  the  carbon-Pt 
mixture.  This  can  be  seen  more  clearly  in  Fig.  2(e),  with  the  lighter 
areas  representing  an  area  of  the  carbon-Pt  mixture  that  contains  a 
higher  level  of  Pt,  compared  to  the  darker  areas  of  carbon  containing 
less  Pt. 

The  isotherms  for  the  adsorption  of  nitrogen  on  samples  CU1, 
CL1  and  CHI  are  shown  in  Fig.  3.  The  isotherms  are  all  typical  of 
a  type  II  response,  although  the  amounts  adsorbed  are  very  low. 
Type  II  isotherms  are  usually  the  result  of  monolayer-multilayer 
adsorption  on  the  open  surface  of  a  nonporous  or  macroporous 


Fig.  3.  Nitrogen  adsorption  isotherms  for  carbon  paper  samples  CU1  (A),  CL1  (O) 
and  CHI  (□). 


adsorbent.  However,  due  to  the  low  amounts  adsorbed,  it  can  be 
inferred  that  the  structures  are  relatively  non-porous,  or  show  very 
limited  porosity. 

Sample  CL1,  which  has  a  lower  loading  of  Ptcm-2,  adsorbed  a 
greater  amount  than  CHI.  This  implies  that  the  surface  roughness 
(and  hence  surface  area)  is  greater  with  a  lower  loading  of  Pt  when 
these  two  samples  are  compared.  Fits  of  the  data  to  the  fractal  FHH 
model  [11  ]  and  fractal  BET  model  [30]  are  shown  in  Figs.  4  and  5. 

Table  1  illustrates  the  parameters  obtained  from  the  fits  of  the 
data  of  samples  CU1,  CL1  and  CHI  to  the  BET  [25],  fractal  FHH  [11] 
and  fractal  BET  [30]  models.  It  can  immediately  be  seen  that  the 
surface  areas  obtained  are  all  low,  with  sample  CL1  showing  the 
highest.  The  fractal  dimensions  obtained  from  both  the  fractal  FHH 
[11]  and  fractal  BET  [30]  model  follow  a  similar  trend.  The  fractal 
dimension  for  CU1  is  higher  than  that  for  CL1  and  CHI.  The  fractal 


Fig.  4.  Examples  of  typical  fits  (solid  lines)  to  the  fractal  FHH  [  11  ]  model  to  nitrogen 
adsorption  isotherms  for  carbon  paper  samples  CU1  (A),  CL1  (O)  and  CHI  (□). 
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Ln  [1  =  (P/P0) ] 


Fig.  5.  Examples  of  typical  fits  (solid  lines)  to  the  fractal  BET  [30]  model  to  nitrogen 
adsorption  isotherms  for  carbon  paper  samples  CU1  (A),  CL1  (O)  and  CHI  (□), 

Table  1 

Parameters  obtained  from  fits  of  the  BET  [25],  fractal  FHH  [11]  and  fractal  BET  [30] 
models  to  nitrogen  adsorption  data  from  carbon  paper  samples  CU1,  CL1,  and  CHI 


dimension  for  CL1  is  higher  than  CHI,  implying  that  the  surface  has 
a  higher  degree  of  roughness  over  the  length-scales  indicated.  This 
would  support  the  fact  that  the  BET  surface  area  is  higher  for  CL1 
than  CHI,  and  infers  that  a  lower  loading  of  the  Pt  catalyst  results 
in  a  greater  surface  area,  which  has  implications  on  the  catalytic 
activity  of  the  electrode  during  its  use  [37], 

The  two-component  BET  model,  which  has  been  used  elsewhere 
[32],  was  applied  to  samples  CU1,  CL1  and  CHI.  The  results  are 
shown  in  Table  2  and  illustrate  that  the  model  is  valid  for  all  sam¬ 
ples.  For  sample  CU1,  which  is  predominantly  made  of  carbon  fibres 
with  a  PTFE  binder  the  two  component  constants  Ci  and  C2  are  13 
and  50,  respectively.  The  determined  BET  constant  from  the  original 
BET  model  was  16.  The  results  imply  that  the  strength  of  interac¬ 
tion  between  the  carbon  and  nitrogen  is  more  predominant  than 
that  with  the  binder,  as  the  determined  BET  constant  is  closer  to  Ci . 
The  same  can  be  said  for  samples  CL1  and  CHI,  with  components 
C\  and  C2  being  the  carbon  substrate  and  Pt  catalyst,  respectively. 
The  results  show  for  both  materials  that  C2  is  the  more  predomi¬ 
nant  component.  It  can  also  be  seen  from  Table  2  that  the  fractal 
dimensions  obtained  agree  with  those  seen  from  the  other  models 
applied,  indicating  the  validity  of  the  model  and  accuracy  of  the 
results  obtained. 

Table  2 

The  results  of  applying  a  two-component  fractal  BET  fit,  via  Eq.  (16),  to  microporous 
carbon  paper  samples  CU1  (unloaded,  0  mg  cm-2  Pt),  CL1  0.39  mg  cm-2  Pt)  and  CHI 
(4.7  mg  cm-2  Pt) 

Sample  BET  specific  area  (m2  g 1 )  C  Q  C2  d 

CU1  0.81  ±  0.006  16  13  50  2.76  ±  0.02 

CL1  8.38  ±  0.08  51  13  106  2.59  ±  0.03 

CHI  5.61  ±  0.06  36  13  34  2.38  ±  0.02 


Relative  Pressure,  P/P0 


Fig.  6.  (a)  Nitrogen  adsorption  isotherms  for  particulate  carbon  support  samples 
CPI  (O),  CP2  (□)  and  CP3  (A),  (b)  A  more  detailed  image  of  samples  CP2  (□)  and  CP3 
(A). 


4.2.  Particulate  carbon  support  samples  CPI  -CP2 

Fig.  6(a)  and  (b)  shows  the  isotherms  for  the  adsorption  of  nitro¬ 
gen  on  samples  CPI,  CP2  and  CP3.  It  can  clearly  be  seen  that  there  is 
a  significant  difference  in  the  amount  adsorbed  between  the  three 
samples.  All  samples  can  be  classified  as  isotherm  type  II,  although 
CPI  is  indicative  of  a  sample  with  high  C(BET)  and  a  large  volume  of 
micropores.  Samples  CP2  and  CP3  show  relatively  similar  amounts 
adsorbed,  whereas  sample  CPI  adsorbed  ca.  7  and  3  times  as  much 
as  samples  CP2  and  CP3,  respectively.  The  three  samples  have  been 
analysed  using  the  fractal  FHH  [11]  and  fractal  BET  [30]  models. 
Examples  of  the  fits  are  shown  in  Figs.  7  and  8. 


adsorption  isotherms  for  particulate  carbon  support  samples  CPI  (O).  CP2  (□)  and 
CP3  (A). 
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Table  3 

Parameters  obtained  from  fits  of  BET  [25],  fractal  FHH  [11],  fractal  BET  [30],  fractal  Porod  Law  [38]  and  DFT  [39]  for  particulate  carbon  support  samples  CPI,  CP2,  and  CP3 


5AXS  [38] 

DFT  [39] 


Parameter 


Specific  surface  area  (m2  g-1 ) 


Length  scale  (nm) 
R2 


Length  scale  (nm) 
R2 


Lnq  (nm-1) 
Length  scale  (nm) 
R2 


Cumulative  specific  surface  area  (m2  g-1 ) 


795.7  ±2.2 
204 

2.73  ±0.07 
5.78-48.06 
2.02-16.86 
0.997 

2.78  ±0.03 
7.05  ±46.69 
2.47-16.34 
0.9937 
2.76  ±0.06 
-0.5  to  -2.0 
0.5-2.0 
0.9993 
2.81  ±0.009 


n/a 

n/a 

2.74  ±0.03 
1.03-34.49 
0.37-12.07 
0.9954 


1.90  ±0.01 
-0.5  to  -2.0 
0.5-2.0 
0.9981 
1.90  ±0.01 
125.7 


0.9996 
2.6  ±0.01 
0.89-23.94 
0.31-8.38 
0.9972 
2.56  ±0.03 
0  to  -2.0 
0-2.0 
0.9986 
2.61  ±0.01 


adsorption  isotherms  for  particulate  carbon  support  samples  CPI  (O)  and  CP3  (A). 

Small  angle  X-ray  scattering  data  was  acquired  for  the  three 
samples  and  examples  of  the  fits  to  Eq.  (16)  are  shown  in  Fig.  9. 
Table  3  shows  the  parameters  obtained  from  the  fits  of  the  sample 
data  to  the  standard  BET  model  [25],  fractal  FHH  model  [11],  frac¬ 
tal  BET  model  [30],  and  the  fractal  Porod  Law  [32],  It  can  be  seen 
that  the  BET  surface  area  for  the  three  samples  vary  greatly,  and 
from  the  high  values  of  the  surface  area  and  C  constant  calculated 
for  sample  CPI,  it  can  be  inferred  that  the  sample  is  microp- 


angle  X-ray  scattering  data  for  particulate  carbon  support  samples  CPI  (Q),  CP2  (□) 
and  CP3  (A),  according  to  Eq.  (16). 


orous  (average  pore  size  <2nm).  Sample  CP3  possesses  a  surface 
area  that  is  more  indicative  of  a  mesoporous  solid  (average  pore 
size  2-50  nm).  Table  3  also  shows  that  the  BET  model  was  only 
applied  to  samples  CPI  and  CP3.  The  BET  model  was  not  valid 
for  sample  CP2  and  as  a  result,  the  experimental  isotherm  was 
fitted  to  a  distribution  of  model  isotherms  calculated  by  density 
functional  theory  (DFT),  as  described  elsewhere  [39],  The  theory 
assumes  slit  shaped  pores  and  the  fit  of  the  model  to  experimen¬ 
tal  data  is  shown  in  Fig.  10.  The  model  shows  a  very  good  fit  to 
the  experimental  data,  and  implies  that  sample  CP2  is  comprised 
of  micro-,  meso-  and  macropores,  according  to  IUPAC  definitions 
[40], 

The  most  striking  result  from  the  application  of  these  models 
to  the  data  is  that  sample  CP2,  which  exhibited  a  negative  C  con¬ 
stant  from  the  fitting  of  the  original  BET  model,  also  resulted  in 
a  fractal  dimension  obtained  independently  by  SAXS  that  did  not 
match  with  that  obtained  from  the  fractal  FHH  model.  The  other 
two  samples,  CPI  and  CP3,  show  extremely  good  agreement  from 
the  use  of  all  the  models,  especially  CP3  which  resulted  in  a  frac¬ 
tal  dimension  that  was  almost  identical  for  all  the  models  used. 
This  result  strongly  illustrates  the  importance  of  applying  various 
structural  models/techniques  to  investigate  the  characteristics  of 
complex  carbon  materials. 

In  addition,  it  has  been  previously  shown  [32]  that  the  presence 
of  adsorbed  chemical  species  on  the  surface  of  porous  catalysts  may 
have  a  significant  effect  on  the  catalysts  adsorption  properties.  This 
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Fig.  10.  Nitrogen  adsorption  (+)  and  DFT  (o)  isotherms  for  particulate  carbon  sup¬ 
port  sample  CP2  (Vulvan  XC72R). 
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was  shown  to  be  the  case  during  the  adsorption  of  nitrogen  on  a 
variety  of  porous  silica  catalyst  pellets.  It  is  possible  that  a  similar 
effect  occurs  within  sample  CP2,  and,  as  a  result,  an  incorrect  fractal 
dimension  is  obtained. 

5.  Conclusions 

(1)  The  results  obtained  have  shown  that  the  surface  structures 
of  the  three  carbon  paper  samples  (CU1,  CL1  and  CHI)  vary 
significantly.  Sample  CU1  showed  a  low  porosity,  low  surface 
area  and  is  highly  disordered,  whereas  CL1  had  a  high  sur¬ 
face  area  (in  relation  to  the  samples  examined  here)  and  a 
more  increased  order,  whilst  CHI  showed  a  lower  surface  area 
and  fractal  dimension  than  CL1.  In  terms  of  BET  surface  area, 
it  would  be  expected  that  the  samples  would  increase  in  the 
order  CU1  <  CL1  <  CHI.  However,  this  was  not  the  case,  perhaps 
because  of  the  large  cracks  present  in  CL1,  which  would  result  in 
an  increased  adsorbance  of  adsorbate.  In  terms  of  fractal  dimen¬ 
sion,  the  results  obtained  followed  an  identical  trend  to  that  of 
the  BET  area.  If  the  cracks  seen  in  CL1  were  not  present,  the 
BET  area  would  decrease,  and  the  resultant  fractal  dimension 
would  be  less  than  that  of  CHI.  The  cracks  seen  in  CL1  were 
formed  due  to  stress  experienced  through  the  Pt  layer  during 
the  deposition  process. 

(2)  Three  porous,  particulate  carbon  blacks  (CPI,  CP2  and  CP3) 
used  as  electrocatalyst  supports  for  MEAs  have  been  charac¬ 
terised  using  gas  adsorption  and  a  variety  of  structural  models. 
The  results  have  shown  that  sample  CPI  was  highly  porous, 
whilst  sample  CP3  was  relatively  mesoporous.  The  fits  of  both 
samples  to  the  fractal  FHH  and  BET  models  showed  very  close 
fractal  dimensions  and  implied  that  the  models  are  valid  for 
the  materials.  The  fractal  dimensions  obtained  from  the  two 
models  agreed  within  experimental  error  to  the  fractal  dimen¬ 
sion  obtained  from  SAXS.  These  results  implied  that  the  models 
accurately  describe  the  surface  roughness.  Sample  CP2,  how¬ 
ever,  did  not  show  a  valid  C  constant  derived  from  the  BET  model 
and,  as  a  result,  the  fractal  BET  model  was  not  valid.  The  sur¬ 
face  area  values  calculated  by  DFT  imply  that  the  material  was 
comprised  of  micro-,  meso-  and  macropores. 

(3)  The  models  used  were  valid  for  samples  CPI  and  CP3  but  not 
CP2  The  behaviour  of  CP2  may  be  explained  by  the  presence 
of  chemical  species  on  the  pore  surface,  which  may  interact 
with  the  adsorptive  leading  to  a  complex  surface  layer.  The  fact 
that  the  BET  model  did  not  hold  for  CP2  may  indicate  certain 
restraints  in  the  assumptions  of  the  model,  namely  that  the  BET 
model  assumes  chemically  heterogeneous  surfaces.  The  results 
for  samples  CPI  and  CP3  suggest  that  the  models  used  were 
valid,  as  the  results  obtained  from  gas  adsorption  agree  with 
those  obtained  from  SAXS.  This  is  apparent  from  two  models 
with  different  theories  applied  to  two  different  physical  pro¬ 
cesses. 

(4)  The  results  obtained  here  can  be  used  to  provide  information 
on  the  design  of  future  PEM  electrodes  and  MEA  assemblies.  By 
understanding  the  physical  characteristics  of  existing  materi¬ 
als,  such  as  those  examined  here,  in  more  detail,  it  is  possible 


to  fabricate  future  materials  that  may  improve  the  operating 
performance  of  PEM  fuel  cells. 
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